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Abstract 

In the theory of random matrices, several properties are known to be uni- 
versal, i.e. independent of the specific probability distribution. For instance 
Dyson's short-distance universality of the correlation functions implies the 
universality of P{s), the level-spacing distribution. We first briefiy review 
how this property is understood for unitary invariant ensembles and consider 
next a Hamiltonian H = Hq + V , va. which Hq is a given, non-random, N by 
N matrix, and V is an Hermitian random matrix with a Gaussian probabil- 
ity distribution. The standard techniques, based on orthogonal polynomials, 
which are the key for the understanding of the Hq = case, are no longer 
available. Using then a completely different approach, we derive closed ex- 
pressions for the n-point correlation functions, which are exact for finite N. 
Remarkably enough the result may still be expressed as a determinant of an 
n by n matrix, whose elements are given by a kernel K{\, fi) as in the Hq = 
case. From this representation we can show that Dyson's short-distance uni- 
versality still holds. We then conclude that P{s) is independent of Hq. 
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I. INTRODUCTION 

Many years ago Wigner [|I| introduced the level-spacing probability distribution P{s), in 
his discussion of nuclear energy levels. The exact form of P{s) was found later in the theory 
of random matrices for the Gaussian unitary ensemble (GUE) ||^-§]. This level-spacing 
probability distribution P{s) was empirically found to be universal in many different cases, 
for instance non-Gaussian probability distributions, or band matrices (in which case the 
measure is not unitary invariant), and even for problems of an a priori different nature 
such as the level spacing of the zeros of the Riemann zeta function 0,^,^which is known to 
coincide with that of the GUE. 

In the next section, we first review how the universality of P{s) has been derived for 
non-Gaussian unitary invariant ensembles, in which the probability measure is given by 

Zj 

where /(x) is an arbitrary polynomial. One first integrates out the unitary group, in order 
to obtain a probability distribution for the eigenvalues of H . It is then easy to show that 
the n-point function may be written as an n by n determinant; the matrix elements of this 
determinant are given by a kernel expressed in terms of orthogonal polynomials with respect 
to the weight exp[— iV/(x)]. Then the understanding of the relevant asymptotic behavior 
of these polynomials at large order allows one to prove the short-distance universality of 
this kernel. From there one can derive the universality of -P(s) in the scaling limit in which 
A^ goes to infinity, the distance x between two neighboring eigenvalues, goes to zero, and 
s = Nx is held fixed. 

In the third section we consider a Hamiltonian which is the sum of a given determin- 
istic part Hq and of a random potential V with a Gaussian probability distribution. The 
measure is not unitary invariant, but one can still write the probability distribution for the 
eigenvalues of H through the well-known Itzykson-Zuber integral |^. Generalizing a method 
introduced by Kazakov ||^ for the density of eigenvalues, we write a representation of the n- 
level correlation function, in terms of an exact and explicit integral over 2n variables. Then 
one discovers that an amazing algebraic structure allows one to express again this n-point 
function in terms of a determinant of an ra by ra matrix. The matrix elements are given by a 
kernel which has an explicit representation as an integral over two variables. In a previous 
paper 0, we had discussed already the two-level correlation function of this Hamiltonian 
through the same method, and we had shown that the behavior of this correlation function 
is indeed universal, i.e. independent of the Hamiltonian Hq, in the short range scaling limit, 
in which the distance x of the two energy levels becomes small, and N goes to infinity, with 



fixed Nx. We had also briefly discussed the n-point function in |TD|. The main steps are 
recalled here; the universality of P{s) follows immediately. 

In the last section we establish some properties of this kernel, and show that it does 
satisfy some necessary consistency conditions. 



II. LEVEL-SPACING DISTRIBUTION P{S) FOR GENERALIZED GUE 

ENSEMBLES 

We return to the single random matrix case with a probabihty 

P^H) = le-^^^(^) (2.1) 

Zj 

and integrate out the unitary degrees of freedom. The resuhing probabihty distribution for 
the N eigenvalues of if is 

PAr(xi, . . . , xat) = C n(x. - x,)2e-^St, /(-») (2.2) 

The n-point correlation function Rn{xi, ...,a;„), is defined as 

JYJ roo /-oo 

i?„(a;i,...,a;„) = -— — '-^ ■■■/ dxn+i- ■ ■ dxNPN{xi,- ■ ■ ,xn) (2.3) 

(A/ — raj! J-oo J-oo 

Following Mehta 0, we introduce the orthogonal polynomials (j)k{x) with respect to the 
measure exp[— A^/(a;)]. Then i?„ is given by the determinant, 

Rn{Xi, ■■■,Xn) = det[KN{Xi, Xj)]ij=i,...,n. (2.4) 

in which the kernel K]\f{x,y) is expressed as a sum of orthogonal polynomials 

1 N-l 

KNix, y) = e-f (/(-)+/(^)) ^ Ux)Uy)- (2.5) 

^^ fc=0 

For instance the pair correlation function, the n = 2 case, becomes 

R2ixi,X2) = pixi)p{x2) - Kn{xi,X2)Kn{x2, Xi) (2.6) 

in which the density of states p{x) is the diagonal part of the kernel p{x) = K]\f{x, x). With 
our normalization conventions the density of state p{x) has a support of finite extension in 
the large N limit. 

In the short distance scaling limit, Kn{xi,X2) becomes Pfl,]!]]] 

^ , , sin[vriV(xi - x,)p(fel±^)] 

^"("^'"^)" nNix,-X2) ^'-'^ 



for A^ -^ oo, xi — a;2 -^ and finite N(xi — X2). The universality of ( p.TD with respect to the 
function /(x) which characterizes the probability measure is thus manifest. The universality 
of the level-spacing distribution P{s) follows at once. 

Indeed, following Mehta @], we first compute the probability E{6) that the interval 
[— |, |] does not contain any of the points xi, ...,xn- in the large N limit. It is thus obtained 
by integrating the N variables of Pn{xi, ...,xn) outside the interval [— |, |] : 

P'i^) = / ■■■/ Pn{xi, ...,XN)dxi- ■ ■ dxN (2.8) 

J out J out 



where the integrals are performed outside the region [— |, §]; 

f dx={r - f^)dx (2.9) 

J out J — CO J — -7r 



We may thus express E{6) in terms of the RnS by using systematically ( |2.9| ) for all the N 
variables: 



2 r<L ,.1 



E{e) = l-Nj_^ p{x)dx + ^J_J_, ^2(x, y)dxdy + ■■■. (2.10) 

The natural scale for the level spacing 9 is of order -^ since in the large N limit the support 
of the density of state is finite. We thus consider the short distance scaling limit, in which 
6 goes to zero and N to infinity, with fixed A^^. In that scaling limit 



N / ^ p{x)dx = NOpiO) + 0{1/N). (2.11) 

We thus define the scaling variable 

NOpiO) = s. (2.12) 

The next terms of ( p^.lO| ) are obtained in this limit by the change of variables 

Nxp{0)=x'. (2.13) 

Then, in the scaling limit. 



s „ s 



N^J] J] R2{x,y)dxdy = j\ j\ ^x' ,y')dx'dy' (2.14) 



2 2 2 2 



with 



in which 



Rn{Xi,---,Xn) =det[k{Xi,Xj)]ij=i^,„^n- (2.15) 



f>f N smvr 7/1-7/2 

K{yi,y2} = 7 — . (2.16) 

7r(l/i - 1/2) 

In this scaling limit we thus obtain 

E{s) = J2 i~ / ' ■■■ /' dxi, ■ ■ ■ , dXndet[K{xi, Xj)]ij=i^,„^n (2.17) 



n=0 

From this representation it is easy to expand E{s) for s small; for instance 

36* 675 



II II ^2(x,i/)rfxrfy = ^s' - ^s' + 0(.«) (2.18) 
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and since the n = 3 term of (|2.17|) is easily shown to be of order s^ for s small, we find 



2 4 

E(s) = l-s + — s^-— s^ + Or/). (2.19) 

^ ' 36 675 ^ ^ ^ ' 

One can also introduce the eigenvalues Aj(s) of the integral equation for the kernel K on 
the interval [— f , +f], 



_^ K{x,y)^i{y)dy = XM^)- (2-20) 

2 

From (|2.17| ) we can write 

oo 

Eis) = n(l - A.) = det[l - K]. (2.21) 



i=l 



For small s, a perturbational expansion using Legendre polynomial gives the same result as 
( p.l9| ). The level-spacing probability distribution P{s) is now obtained from E{s) 

P{s) = ^^E{s) (2.22) 

Through this representation, we find that the universality of P{s) results from two 
sources; i) the n-point correlation Rn is expressed as the determinant of a kernel K]^{xi, Xj), 
ii) the kernel Kn{x, y) has a universal short distance behavior K in the short distance scaling 
limit. 



III. DETERMINISTIC PLUS RANDOM HAMILTONIAN 

We now consider a Hamiltonian H = Hq + V , where Hq is a given, non-random, N x N 
hermitian matrix, and V^ is a random Gaussian hermitian matrix. The probability distribu- 
tion P{H) is thus given by 

P{H) = ie-f->'' 

= _Le-f i««»-2«o«) (3.1) 

Zj 

We are thus dealing with a Gaussian unitary ensemble modified by the external matrix 
source Hq^ which breaks the unitary invariance of the measure. In previous work P,pD|, 



we have already discussed the density of state, and the two-level correlation function. For 
completeness, we repeat here the basic steps. The density of state p(A) is 

p[X) = 1 < Tr<5(A -H)> 

+ 00 rU 

— e-*^*"f/(t) (3.2) 

oo 27r 



where U(t) is the average "evolution" operator 

U{t) =< Tie'''''' > . (3.3) 



We first integrate over the unitary matrix u! which diagonahzes H in ( p.l| ), and without 
loss of generality we may assume that Hq is a diagonal matrix with eigenvalues (ei, ■ ■ ■ , cn)- 
This is done with the help of the well-known Itzykson-Zuber integral |0], 

/..e.p(T...B..,^M!^ (3.4, 

where A{A) is the Van der Monde determinant constructed with the eigenvalues of A: 

N 

AiA) = l[ia,-a^). (3.5) 

i<j 

We are then led to 

1 ^ r 

^(^) = y'Afu^ E / ^^1 ■ ■ ■ dzNe'^^'^-Aixi, ■ ■ ■ , x^) 

xexp(-^^x2 + iV5]e,x,). (3.6) 

The normalization is fixed by 

t/(0) = A^ (3.7) 

The integration over the Xj's may be done easily, if we note that 

f N 

dxi-- ■dxNA{xi,-- ■,XAr)exp( — 2^^1 + N^hiXi) 

= A(6i,---,Mexp(yE&') (3-8) 

Putting hi = ei + it6a,i, we obtain 

f/(t) = E n ( '2 '^ )^-^+mte^ (39) 

The sum over N terms in (|3.9| ) may then be replaced by a contour-integral in the complex 
plane, 



1 f du -r-r ,u — e^ + it 



2 



^w-s/iini^ir^)--^-'"" (3,10) 

The contour of integration encloses all the eigenvalues e^. The Fourier transform with respect 
to t gives the density of state in the presence of an arbitrary external source Hq and for 
finite A^. 
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In the case of the two-point correlation function, we have 

i?2(A, /.) =< ^Tr(A - H)^TTif, - H) > (3.11) 

By using integral representations for the two 5-functions, the two-point correlation function 
i?2(A,/i) is expressed as the Fourier transform of U{ti,t2), 

UitiM) =< Tre*^*^^Tre*^*^^ > (3.12) 

Using again the Itzykson-Zuber formula to integrate over the unitary group, we obtain 

TV „ N /\(t\ 



After integration over the Xj's, we have 



Y\i<j{(^i - gj + ^^l(^»,ai - ^i,ai) + ih{^i,a2 ~ ^^,^2)) 



AT ,2 N ,2 



^ Afitieai+Afit2<:Q2~"T*l~T*2~''^*l*2'5ai,Q2 ('3 \A) 

This term is devided into two parts; ai = 0:2 and ai 7^ 0^2 cases, 

f/(tl, t2) = yn ^^' ~ ^2 + ^(tl + t2)(^»,ai - ^i,aj ^jv^fa+t2)e^,-^fa+t2)^ 

(e^i - Caj -hi(ti - ^2)) -p-r (cai - e^ + iti) (caj - e^ -^2^2) 



^~^ \^ai "-a2 ' "V"! ''^// TT 



.iV/.2 I i2 



^gAfJtl<:Qi+Af*t2<:Q2-f (*?+*!) ("3 ^5) 



Fourier transform of the first term becomes 5-function ||T4[, and can be neglected for -R2(A, /x) 
for A 7^ /i. The double sum in ( |3.15 ) may be written again as an integral over two complex 
variables: 

U(t^ t^) = -J—e'^'l'^'l / dudv ^^a^^^N^t2v {u-v+ {ih - it2))iu - v) 
(tit2) J {27riy {u — V + iti){u — V — it2) 



^ ih . , iU 



X 



Noting that 



^=1 {u - e.f) {V - e^ 



^ ^1^2 ^ {u-V + i{ti -t2)){u-v) 

{u — V + iti){u — V — it2) {u — V + iti){u — v — it2) 



we find that (p.l6|) is a sum of the disconnected term and a connected part. We know 



Fourier transform U with repect to ti and ^2 and shift the integrations variables. By the 
shifts ti — > ti + iu, and ^2 -^ t2 + iv, we see easily that -R2(A, /i) is a two by two determinatnt, 
namely that 



i?2(A, /i) = Kn{X, X)KN{fi, /i) - Kn{X, iJi)KN{ii, A) (3.18) 

with the kernel 

^^(A,/.) = /^ /i^ n(^^^^)^e-f^^-f*^--— '^ (3.19) 

J In J iTti ^\ t> — e^ V — it 

Note the similarity of the determinantal structure found here with that of the zero source 
case given in ( p.4| ). 

In [^, this kernel Kn{X, /j,) was examined in the scaling limit, large N, but fixed N{X — fi). 
In this limit one can evaluate the kernel ( p.l9| ) by the saddle-point method. The result was 
found to be, up to a phase factor that we omit here, 

Kn{X^, A2) = -— sin[7ryp(Ai)] (3.20) 

Try 

where y = N{Xi — A2). Apart from the scale dependence provided by the density of state 
p, the two-point correlation function has a universal scaling limit, i.e. indepent of the 
deterministic part Hq of the random Hamiltonian. 



IV. DETERMINANT FOR THE N-POINT CORRELATION FUNCTION 

The n-point correlation function -Rn(Ai, ■ ■ ■ , An) is given by 

i?„(Ai, ■ ■ ■ , A.) = -i^ < n Tr5(A, - M) > (4.1) 

If we put the constraints that all Aj are different, this expression conincides with (|2.1| ). When 
some Aj become same, we have extra 5- functions as shown in ^. Therefore, we assume all 
Aj are different. 

Without an external source, this n-point correlation function is expressed in terms of the 
kernel K^iXi, Xj) as 0,|[ 

i?n(Ai,---,A„) = det[ir^(A„A,)] (4.2) 

where i,j = 1, ■ ■ ■ , n. This result was derived by the use of the orthogonal polynomials. In 
the external source problem, we can not apply the orthogonal polynomial method. Our aim 
is to find a proof of ( [4.2|) for the external source case. 
Using the Itzykson-Zuber formula of ( p.5|) , we have 

i? (Ai ■■■ A ) = — V / ^^i---^^»^(-g) gfE''?+»E*^A, .4 3^ 



Oij^a 



where 

bk = ek + i{tlSk,a^ H h tnSk,an)- (4-4) 



Using the contour-integration representation, we get 

" ~ ' (27r)" 

dui---dun myt^u, rr TT n + 



X 



itr 



^ TT 1 -pr K - ^g + ijtp - tg)]{Up - Uq) 
p P p<q \^P ^g "•" ^'^p)\^p ^q ^''q) 



(4.5) 



When n = 2, this reduces to the previous expression p.l6| ). We make a shift of the variables 
tp-. tp -^ tp + iup then we get 



Rn 



(27r)" / {2'KiY 



n N 

nn 



. —ta + it 



a I "''p ^ 



n( 



Xlp %bq , [Up Uq) TH 



n 



We recognize in ( [4 .61 ) a Cauchy determinant, 
det 



1 



lOj - bj 



n(n-i) n^<j(ai - aj)(6i - 6j) 



Uu{ai-bi) 



i,j=l,...,n -^-^^J* 

if we identify ak to itfe, and bk to Wfc in (|4.6| ). Then, Rn is given by 



Rn 



(27r)" 7 (27ri)" 






fc=i «=i ^o^ ^fc * ^i 

(27r)" 7 (27ri)" 



X det 



AT 

n 



'iZi + €qi 



Using the expression for the kernel of ( p.l9|) , we obtain 

i?„(Ai, ■ ■ ■ , A„) = det [Kn{\u -^i)]ij=i,„,,n 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



We could thus prove the determinantal form of the n-point correlation function for a deter- 
ministic plus random Hamiltonian. 



V. THE PROPERTIES OF THE KERNEL 



As we have seen in ([4.9| ), the n-point correlation function Rn is expressed by the determi- 
nant in the presence of the external source. When we integrate out the variables a;/+i, ..., a;„ 
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of i?„(a:i, ..., x„), we obtain the 1-point correlation function Ri{xi, ...,xi). Since we have ([4. 91), 
the necessary consistency condition for this result is 

dfiKN{X, fi)KN{fi, u) = K^iX, y) (5.1) 

This property is verified easily by the contour-integral representation of the kernel i^Ar(A, /x) 
given in ( glSj) 0. We have 



Kn{\, h)Kn{ijI', i')dfj, 



oo 
oo 



dtidt2 [ dUidU2 YT/^l + ^'tl^,^'y + ^2^ 1 

(27r)2 / (27ri)2 -1-^^^ ^^ _ e^ ^ ^ „2 - e^ (wi + «ti)(M2 + ^ta) 

^^-^;{ul+ul+tl+tl)-^Nt^\~^Nt2^Ji~Nu^^JL■-Nu2U /^ 2) 

Integration over yU, after the shift ^2 ^ ^2 + iwi, gives a delta-function for t2, and the 
contour-integral over ui around the pole ui = —iti reconstructs Ki^{\, v). 

We observe also the kernel KN{X,fi) has A^ eigenvalues equal to one, with Hermite 
polynomials as eigenf unctions since, for n < N, 

Kn{\ ^)Hn{\fNii)t--^ dfx = HniVNX)e--^ (5.3) 

-oo 

with, Hq{x) = 1, Hi{x)= X, H2{x) = x^ — 1, etc. This property may also be easily verified 
through the contour integratal representation. For n > N — 1, ( p. 3D does not hold. The 
right hand side of (|573|) becomes e^ dependent. When the external source e^ goes to zero, 
the right hand side of ( |5.3| ) is vanishing for n > A^ — 1. This is of course related to the fact 
that the kernel is then expressed as a finite sum of Hermite polynomials. 



VI. SUMMARY AND DISCUSSION 

In the previous section, we have proved that the n-point correlation function is expressed 
by the kernel K{x, y) as in the absence of an external source. In the short distance limit, in 
which (A — ^)N is kept fixed, the kernel Kisi{Xp, Xg) takes a universal form, and the n-point 
correlation becomes universal ( up to a rescaling by the density of state p). As we have seen 
in section 2, the level-spacing probability distribution P{s) is given by an integration over 
the n-point correlation function Rn{Xi, ..., A„). Therefore, in the short distance scaling limit, 
P{s) has a universal form, independent of the deterministic part. 

We have assumed throughout this work that the density of state is finite,of order one, in 
the energy range that we are considering. We have discussed the level-spacing probability 
P{s) for two levels centered around the enregy zero. If we considered instead two levels 
centered around an energy Eq, i.e. an interval [— | + Eq, | + Eq], the behavior of the kernel 
K]\r{X,fi) remains universal, apart from the scaling by the the density of state p(ii^o)instead 
of p(0). Therefore, we still have the same universal spacing dsitribution P{s) for an arbitrary 
energy Eq as long as p{Eq) remains of order one. 
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We have also assumed that the eigenvalues of the deterministic term Hq are inside the 
support of the asymptotic smooth density of state p. When the eigenvalues are widely 
separated, the density of state shows an oscillatory behavior. In such cases, the two-level 
correlation function, or the kernel K]y{X, /x) does not approach, in the scaling limit, the sine- 
kernel, and a universal form for P{s) is not expected |T^. This is reasonable, since we know 



that, when the random potential V increases in comparison with the unperturbed determin- 
istic term Hq, we crossover to a universal behavior independent of the initial deterministic 
term. 

Finally we have found two kinds of universality: either Hq = and the distribution of 
H is non-Gaussian, or Hq is non-zero and V is Gaussian. It is tempting to conjecture that 
this generalizes to non-Gaussian problems with a non-zero source as considered in the case 



of the two-level correlation function [|1T]-IT^ , or to the time-dependent case [IIO 
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